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1 Introduction 



In recent years, relativists have been researching on what happens when 
gravity is coupled to nonlinear field theories, as for example the Yang-Mills 
(YM) fields. In the beginning, there was doubt about whether there were 
any solutions, since neither the vacuum-Einstein equations [[jj, nor the pure 
YM equations have nontrivial, static, globally regular solutions [Q, [|. 

Thus, it was even more surprising, as in 1988 Bartnik and McKinnon (BK) 
discovered soliton solutions to the SU(2) Einstein- Yang-Mills (EYM) system 
numerically [|J. These solutions provided for the first time solid evidence 
that the EYM system has particle-like solutions. Later, other authors M, 
H [7[] found the colored black hole solutions to the same model, which are 
a counterexample to the no-hair conjecture; i.e., the well-known uniqueness 
theorem for stationary black holes of the Einstein-Maxwell system has no 
natural generalisation on the non-Abelian case. Soon, the existence of both 
types of solution was established rigorously |8|, |, [10], |TTJ, and their instability 
was proven |12|, [13|, [U], [15| . Because of this instability, it was expected that 
it would also apply to arbitrary gauge groups. Proof was first successfully 
provided for the EYM solitons |16| and shortly afterwards for black holes as 
well |TjJ. 

There were also attempts to find corresponding solutions for other re- 
lated field theories. The Einstein-Skyrme system for example has black hole 



solutions with hair that are at least linearly stable |1| H91 GO, 21]. Several 



authors looked at other models, notably the SU(2) Einstein- Yang-Mills-Higgs 
(EYMH) theory with a Higgs triplet ^ as well as the EYM-dilaton 

theory [p5[| . Recently, the instability of the gravitating, regular sphaleron 
solutions to the SU(2) EYMH system with a SU(2) Higgs doublet, that were 
known numerically for some time [pTfl , could be shown [f26|j . 

In November, 1993, Ding and Hosoya (DH) presented a spherically sym- 
metric, analytic solution to the SU(2) EYM system with a positive cosmo- 
logical constant (EYMa) [p8[| . Although, it turned out that the DH solution 
is not really new, but gauge equivalent [^9] to the cosmological solution of 



[ pOfl , it triggered the work presented in this paper. It is natural to expect 
that the cosmological solution given in |30[ is just the ground state of a dis- 
crete family of static spherically symmetric solutions living on a topological 
3-sphere. For reasons to be explained later, one will, however, have to adjust 
the cosmological constant for each member of the family. In my paper I will 
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present the first excited solution which I found numerically. As in the BK 
case, further solutions were expected. In contrast to the analytic solution of 
Hosotani or Ding and Hosoya, g 00 is no longer constant in my solution. Fur- 
thermore, I will prove some general properties of the regular solutions, e.g. 
the interesting fact that solutions on S* 3 can only exist for a non-vanishing 
cosmo logical constant. 

The paper will be organized as follows: In chapter 2, I will first show a 
derivation of the SU(2) EYMa system and the simplest way for obtaining 
the DH solution. This will be followed by the above-mentioned general pro- 
perties of regular solutions. In chapter 3, I will tell how I found the solutions 
numerically, while the conclusions and the discussion can be found in chap- 
ter 4. In conclusion, I will present a cosmological embedding of the EYMa 
system which remarkably simplifies the consideration of stability of the DH 
solution. The idea can be attributed to N. Straumann. 
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2 The EYMa differential equations 

2.1 Derivation of the differential equation system 

We write the metric for a static, spherically symmetric spacetime in the usual 
way: 

g = - e 2a dt 2 + e 2b dr 2 + r 2 (d$ 2 + sin 2 ?? dip 2 ) . (1) 

a and b are functions depending exclusively upon r. 

We are only interested in YM fields that are purely magnetic. According 
to Ref. |17|1 ) we can determine the SU(2) gauge potential A as 

A = (—w T\ + w r 2 ) d$ + (r 3 cot $ — w T\ — w r 2 ) sin# dip . (2) 

w, w are dependent upon r, and r 1; r 2 , r 3 are the spheric SU(2) generators 
with f — a/2i (B being Pauli matrices). With a gauge transformation we 
can furthermore make w disappear, with A finally assuming the following 
form: 

A = w t 2 di) + (r 3 cos $ — w T\ sin ■&) dip . (3) 

This means we are looking for a = a(r), b = b(r) and w = w(r). 

With (0), we can calculate the components of the field strength tensor 
F = dA + A A A 



F 12 = w'r 2 e~ b - , F 13 = -w'n e- b - , F 23 = (w 2 - 1) r 3 \ . (4) 



We made use of the often employed orthonormal basis of 1-forms 

6° = e a dt , 6 1 = e b dr , 6 2 = r dd , 6 3 = r sintf dip , (5) 
and the notation '= d r . The definition of the YM-Lagrangian leads us to 

L YM = -^tr(F^), 



L - = - W ~^2~ + TKZ ,4 • (6) 



1 /2 e- 2b 1 (w 2 -l) 2 

57r r" 1 lo7r r 4 

In the General Theory of Relativity (GR) we obtain the Einstein field 
equations by varying the action 



S 
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with respect to the metric: 



rpflU 

Atx 



F F a 

1 an 1 v 



(8) 
(9) 



The components of the energy-momentum tensor become 



Too 
T n 
T 22 = T 33 



26 

87T ^ r 2 16-7T 
1 (w 2 --!) 2 

16tt r 3 



16tt r 4 
1 ( W 2 -l) 2 



(10) 



and the G^ relative to the basis 



are 



u o 
G\ 

(_r 2 — Lr 3 



+ e" 2b (^ - f ) , 
^2 "I - e (^2 H jr) > 
e~ 26 (a' 2 - a'6' + a" + ^) 



(11) 



All other equal zero. 

The spherically symmetric Einstein field equations with A-term are 



We insert (0) and (0) into (0): 



-26 



1 -2b f 1 2fl \ , , 

- — + e 2b (— + — ) = Gw 



,2 e 



-2b 



G (w 2 - l) 5 



G {w 2 - If 



A 



- A 



(12) 



(13) 



(14) 



For the YM equations, we can operate in a similar way as in obtaining in 
the GR the field equations from a variational principle. Here, the action is 

S ym=J^ymV ■ (15) 

If we now vary with respect to the gauge potential A, we obtain 

D*F = . (16) 
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In addition, we can also prove the identity 

DF = D(dA + A A A) = 0. 



(17) 



The best way to do this, is to use D . = d . + [A, . ] , as well as 
[B, C] = B A C - (-l) Pq C AB for Be A P (M), C G A«(M). The equa- 
tions ([TJ]) and (|I7D are the known YM equations. 

But now, back to our problem. We write the equation for our A 
and our basis (|5|) in full: 



-2b 



w" + w'(a' - V) (w 2 - l)w 



. 



(18) 



(|18D is the equation of motion of w. The equations fll3D , and (|T^) 
constitute a complete set of equations for the static spherically symmetric 
SU(2) EYMa system. Now, we have three differential equations for a(r), 
b(r) and w(r). We write them down for m(r) and S(r) instead of for a(r) 
and b(r), with 



S(r) 

r 2b(r) 

N(r) 
S(r) 



—a(r) — b(r) , 
2 m(r) 

1 



-2b(r) 



-S(r) 



1 - 



2 m(r) 



From (0)-(0) and 5' 



(a' + b'), we obtain 

t 
r 



(19) 
(20) 

(21) 
(22) 

(23) 



Then we differentiate p0| ) with respect to r and insert the result into ( PD 
(with regard to ( ^I|) ) : 



We set I/: 
obtain 



(«. 2 ^-i) 2 

2r 2 ) 



G l2 A _. . G (^ 2 
— ty' iV(r) H — 



A 2 
— r 
2 



(24) 



2 4 
note that N' = ^ - — , multiply both sides by 2 and 



2m' = G(V + Nw ,2 )+Ar 2 
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With 2m! = ^ — rN' and iV — 1 = — — the whole turns into 

r r 

-rN' = G(V+Nw' 2 ) + (N- 1) + Ar 2 . (25) 
We write (|THj) in a slightly different way 

w 2 - 1 

e~ 2b w" + e~ 2b w' (a' - V) - w = 0. (26) 

Then, we consider the expression (^Ne~ s w'^j . Writing it out and multiplying 
by e 5 gives 

(iVe" V) ' = e~ 2b w' (a' - b') + e" 26 w" . (27) 
Comparison of (|27| ) with (p6|) gives us: 

= 2^1 ««-». (28) 

The definition of V leads to: 

1 dV _ w(w 2 - 1) ^ = y 

2 <9w r 2 dw 

Thus, ( p5[ ) becomes: 

(iVe-V)' = iy, w e- 5 . (29) 

We summarize the results: 

— r N' = G {V+Nw' 2 ) + (iV — 1) + Ar 2 (30) 
- r( J' = G/ (31) 



(iVe-V)' = ^e" 5 , (32) 



with 



2 



(1- W 2 ) 2 

1 W(W 2 — 1) 



2 r 

(MD " (H) constitute for the functions N(r), 5(r), w(r) our complete system 
of differential equations of the SU(2) EYMa system. 



2 
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2.2 An analytical solution to the system 



In 1993, S. Ding and A. Hosoya presented an analytical solution to the system 
fl30D - (|32|) [28|. I will give within the framework of our formulas another 
derivation. 

We make the ansatz S(r) = N(r) a , a = —\. This, along with (p2|) , gives 

us 

1 N' 

S '-2N- < 33 » 



Our system (j30|) - (|32|) thus becomes 



From 



-rJV' = N-l + G{V + Nw') + Ar 2 

-rN' = 2GNw' 2 , 
1 



2VN 



v;, 



— (|35|), we obtain 



GNw' 2 = N-l + GV + Ar z 



(34) 
(35) 

(36) 



or 



G(VNw') 2 = N-l + GV + A 



r 2 . 



We calculate the derivative of (|37[) with respect to r: 

2G(v / aV)(v / A%') / = (N-l + GV + Ar 2 )' . 
On the other hand, (|3"5|) leads to 

2G(v / AV)(v / iVV)' = G?V; ifl u;'. 



(37) 



(38) 



(39) 



We compare (|38|) with (p9|) , bearing in mind that V = V lW w' + V^ r 
K,„ w' — - V. Thus, we obtain 



V 



N' = 2G 



In (0), we eliminate i\T' with @ 



A%' 2 + 1/ 



2Ar 



A 
G 



r 2 . 



(40) 



(41) 



We insert this into fl34T) and obtain 



* = -7 + (I ~ 2Ar) • (42) 



We clearly see that ( f4"2| ) is linear. 
The corresponding solution is: 



N= - + (l--Ar 2 
r v 3 

c is an integration constant. With the condition N(0) < 00, c equals zero. 
Thus: 

N = 1 - ^Ar 2 . (43) 
We differentiate fl43|) : AT' = — | Ar, insert this into ([40|) and solve for w: 



3G 

If we insert fl43|) and (fH|) into (|35|) , we see that 



v? = l-\~i Jl . (44) 



1 = 4- ( 45 ) 

According to ([45]), we equate fl4"3] ) with ( PD and obtain a quadratic equation 
in A: 

A(fA-x)-0. 
The sensible solution is A = ^ which results finally in the DH solution 

u, 2 = N = 1 - — . (46) 

Gr 

2.3 Interpretation of the solution 

Our ansatz has been: S = N" 1 ^ 2 , which means that e 2a = 1. With this and 
with fl4Tf ), the metric ([I]) becomes 

3 = -rft 2 + —Lj. rf r 2 + rW . (47) 



We want to replace r by 

u 



1 + — 



and thus re-write our metric: 

l,o u 2 



g = -dt 2 H - g du 2 H 2 dtt 2 . (48) 



Thus, our solution ( fig) describes an Einstein universe [|3T . 
Finally, we observe 

dv 

dg = -r— . (49) 



I 

1 G 



This can be solved: 

r = \/~G sin [ — ^= + c ) . 



c is an integration constant. If we require that r(g = 0) = 0, c becomes zero. 
Thus: 

r = y/G sin-f= . (50) 



'G 

Now, the metric fl47T) can be written as follows: 



g = -dt 2 + dg 2 + G sin2 f-7=J dn2 ■ ( 51 ) 



This means that the spacetime manifold is R x S 3 , g being an angular coor- 
dinate of S 3 . 



2.4 The EYM A equations for L, R, w 

We set up the EYMa equations again, this time for the functions L, R, w of 
the coordinate g. Our metric now reads like this: 

g = -L 2 dt 2 +dg 2 +R\g)dQ 2 , (52) 
R = r , L = VNS . 
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We use our old EYMa equations ( p0|) - (|32| ) and the relation 



We use the notation: 



First, we need 



We insert that into (BO 



dg 



d 
dr 



N' = IR 2 ) = IR 2 ) -^ = 2R, 



d 
dg 

dg 



dr 



w = w 



dg w 
dr ~ R 



Next, we need 



Its derivative results in 



Giy+w 2 ) + (R 2 - 1)+AR' 



S=t = e~ s 
R 



6' 



L R 
~LR R 2 



D and (HD turn @ into 



and solved for L 



With the help of 



LR R 2 R 2 



■ LR „ w 2 L 

L = —. h G r- . 

R RR 
(H) and (H) becomes 

Lw + Lib — — V lW L . 



By inserting ( |59|) into d60|) , we can eliminate L 



.. 1.. w ( „ iir « ^ 
M = _U G hi? 
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And with (|56|), R also disappears: 



1 T . w 
w = - V lW H 

2 ,w 2RR 



G{V-w 2 ) + {R 2 -l)+ARf 



Summarizing: 
R 



2R 



G{V+w 2 ) + {R 2 - 1)+AR? 



w 



w 



2 V ' w + 2RR 
LR „ w 2 L 

+ G r- 

R RR 



G(V-w 2 ) + (R 2 -l) + AR 2 



(62) 

(63) 
(64) 

(65) 



We know from our special solution (^) that R(g) becomes maximal at 
the equator, i.e.: 



R(q) 



. 



To regularize our equations (|63|) - (|65D, we introduce a new function that 
upon considering the equations is perfectly obvious: 



i(q) ■-- 



1 



RR 



G (V- w 2 ) + (R 2 - 1) + AR 



(66) 



Thus, (|63"D, 



and (p5| ) become: 
R = - 



2 ' R 



w = - V rW + - w 7 , 
L = --Lj. 



We also need 7. After a longer, but not at all difficult calculation, we find 
that 



7 = ^7 



R 2 



2GV + R 



S2 



Our final system now looks like this: 



1 1 
w = -V, w + -w-f , 



(67) 
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R 

7 
L 



2 7 ~~ R 2 
2 ' ' 



2G^ + #' 



with 



7 
V 



RR 

(w 2 



G(V-w 2 ) + (R 2 -l)+AR : 
-I) 2 



2R 2 ' 
w(w 2 — 1) 



R 2 



From now on, we will work only with this system. 

In conclusion, I will show how the solution of Ding & Hosoya 
these coordinates. We can find them quickly with the help of 



R 



w 



G sin 



Q 



G, 



cos 



Q 



G, 



L = 
7 = 0, 

q e [o, VGti . 



(68) 
(69) 
(70) 

(71) 
(72) 
(73) 



H looks in 
and (HBfi: 



(74) 
(75) 

(76) 
(77) 



2.5 Expansions 

We already stated that in the case g = we require that R — 0. Furthermore, 
in 7 ([71]) we encounter the term Ai? 2 . As the spacetime manifold has the 
topology of S* 3 , R equals zero again at a certain g. We will call this g . Thus, 
we have 

R(0) = , R(g ) = . 
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For the solution (|74|) - ( fT?|) , g = \fGix. Unfortunately, in our system 
(HI) - flZD ^ a ^ so a PP eares in the denominator; therefore, we have to ex- 
pand our functions R, w, 7, L at g = and g = g - 

Let us first take a look at fl2T|). As we want regular solutions, m(r) — > 
must hold for r — > 0, and thus iV(r) — > 1. According to (|53|) , that means 
that 



# 4 



(e -> or g -> g ) 



UR - 
from 



0, then w — > ±1, since V otherwise becomes singular, which is obvious 
F2|). Equation (|68|) gives us 



RR = RRj - Gw 2 . 

2 ' 

We see immediately that w — > 0, if q — > or q — > g . This corresponds 
to the observation made by BK Q, which says that the solutions become 
asymptotically to w = ±1, since u> must be zero there. 

With enough patience in differentiation and applying the rule of 
Bernoulli-de l'Hopital several times, the equations fl6T|) - (|73| ) provide us 
the important conditions: 



R 



R^O 

-> Gid . 2 

2R 



3R 



g , g -> g . 



R {3> and w can be any values. We set R 3 
We summarize our results again: 



=F2c, c > and w) 



-2b, b > 0. 



i? - 


-> ' 


i? - 


-> ' 


R - 


-> 1 


R - 


-> -1 


R - 


-> 


R - 


-> 


R {3) - 


-> -2c 


>g^0 r {3) - 


-> +2c > 


w - 


-> 1 


w 


-> ±1 


w - 


-> 


w 


-> 


w - 


+ -26 , 


iv 


-> -2b , 



These conditions lead to: 



o 



(78) 
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R = q - ^g 3 , (79) 

w = 1 - bg 2 , (80) 

R = Q?o - e) - | (Qo - q) 3 , (si) 

w = Tl± b(g -g) 2 ; (82) 

± in the formula fl32|) depends upon whether w has an even or an odd number 
of zeros. ( |T8"D - (^) lead to the expansion of 7: 

7 = (-AGtf + ^g, g^O (83) 

7 = (W-jjA) (£?-£?„), f?^^. (84) 
Now we can insert ( |8~3~D and (|84]) into (0), in order to obtain L: 

L = 1 + Q -> (85) 

L = 1 + - ^) (f? - f?) 2 , ^f? . (86) 

As an example, we can expand (ffil) and ( [T5f) and compare them with 

(HID - 0>- Thus, we obtain b = 0.25 and A = 0.75, if we take G = 2. When 
solving numerically, we need those values. 



2.6 Properties of the differential equation system 

We have seen in Sec. 2.5. that the function w starts at 1 and ends at ±1. 
We will now show that w cannot become higher than +1 or lower than — 1 
within the interval ]0, g [, as it would otherwise never return and diverge. 

Furthermore, we will see that any solution with A = is trivial, i.e. 
R = 0, w = ±1. Then I will provide a proof that R can have only maxima 
and no minima within the interval [0, g }. In conclusion, we will recognize 
that at the points where R has a saddle pointQ, w will have an extremum. 
All these facts combined allow for a qualitative idea of the solutions. 

2 point where the first and second derivative vanishes f32| 
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Proposition 1 For any solution w, the following applies: 

Hq)\ < l, Vge]0,g o [. 

Proof: Assume the contrary, i.e., that w is larger than +1 for certain 
f?G]0, g [. In order that the behavior at g = and g = g holds, w must 
have a maximum in that range, where it is larger than +1. That means that 
there is a £>G]0, g [ where the following applies: 

w(g) > 1 , w(g) = 0, w(g)<0. 

Now we use equation (p7|) : 

w (w 2 — 1) 1 

R does not vanish at g and thus, we have 

w(w 2 -l) 

w(g) = 



R 2 
implying 

w(g) > . 

This is a contradiction to the assumption. We can proof the case —1 ana- 
logically. 

q.e.d. 

In order to prove that nontrivial solutions must have A ^ 0, we need several 
tools. 

Proposition 2 Every nontrivial solution has at least one zero of w. Thus, 
there has to be a g* g]0, g [ with w(g*) = 0. 

Proof: We make again a counterassumption and say we have a w for which 
the following holds: 

w(g*)^0, V Q *e]0,Q [. 
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Without loss of generality we can chose w(0) = +1. Then w must have a 
minimum in order to fulfill w(g ) = +1. It cannot be equal to —1, because 
in that case w would have a zero. Therefore, there is a gE ]0, g [ for which 

w(g) > , w(g) < 1 (follows from Prop. P 
w(g) = , w(g) > . 



Taking a look back at equation (|67|) , 

w(g) [w 2 (g) - 1] 



w(g) 



R 2 (q) 



we see that w(g) < 0. 

This is a contradiction to the assumption. We can analogiously treat the 
case u> (0) = —1. 

q.e.d. 

For the ground state with u>(0) = +1, w(g ) = —1, there is exactly one 
zero. 

Proposition 3 // there is a g g]0, g [ with j(g) = 0, then j(g) < 0, if we 
presuppose A = 0. 



Proof: We chose a g with 7(^1) = 0. Equation (|7T|) gives us 
GV(~g) -Gw 2 {g) + R 2 {g) -1 = 0. 



We solve this for G V and insert the result into ( |69D . Thus, we obtain 



i(Q) = -^|y [2Gw 2 {~g] - R 2 {~g) + 1 

Since w 2 and R 2 are always larger than zero and < R 2 < 1 ,7(f?) < 0. 

q.e.d. 

Proposition 4 For every solution with w(0) = +1 and w(g ) = —1 and 
A = 7 must vanish identically. 
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Proof: The equations ([83]) and (|84| ) and the condition A = lead us to 

7(0 = 0) < 

7(0 = 0o) < 0- 

However, since 7(0) = "f(Q ) = still applies, 7 must cross the abscissa at 
least once; e.g. at g, where 7(f?) > 0. This, however, is a contradiction to 
Proposition |3|, which says that 7(0) < 0. Therefore, only 7 = can apply. 

q.e.d. 



Lemma 1 Every solution with A = is trivial, i.e. R = 0, w = ±1. 
Proof: We know from Proposition ^ that 

7 = V e e[o,Q ]. 



Together with ( p9|) and (|7T[) , this leads to 

2GV + R 2 -! = 0, 
GV -Gw 2 + R 2 -1 = 0. 

We subtract the two equations from each other and obtain: 

w 2 = -V . 

Since V > 0, it must be: 

w = 0, V e e[0,Q o ]. 

Therefore, w must be proportional to a constant. In order to keep V regular 
for -R(O) = 0, only w = ±1 can apply. 

With (j68|), we can see that R = 0. Therefore, R should equal zero or be 
proportional to g. However, if it were proportional to g, it could never equal 
zero at g . Therefore: R = 0. 

q.e.d. 

Of course, Lemma [3] only applies if we assume that R equals zero again for 
a finite value of g . But we had already assumed this in Sec. |2~5| . 
Now we will take a look at the properties of R. 
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Proposition 5 The maximum of R is positive, the minimum negative. 
Proof: 

Maximum: 3g* e ]0, g [ , so that 

r(q*) = o, r(q*)<o. 

We insert into equation (|68|): 



R(g*) = -G ' ' 



R{qT) 

w 2 is always positive; therefore, R(g*) > 0. 
Minimum: 3g* e]0, g [ , so that 

r(q*) = o, R{e*)>o. 

With the same argument as above, we show that: R(g*) < 0. 

q.e.d. 

Minima only exist when the radius R at these points is negative. In order 
to fulfill the boundary conditions R(0) = +1, R(g ) = — 1, R must in this 
case have two maxima or two zeros. And for every additional minimum there 
must be an additional maximum. In other words: 

Number of maxima = Number of minima +1. 

If we exclude negative values for all regular solutions to R, we obtain the in- 
formation that all nontrivial, regular solutions to R have just one maximum. 
We also know something about the saddle points. 

Proposition 6 If R has saddle points, then w has an extremum at this same 
point. 

Proof: If R has a saddle point or a flat point, the following applies: 

3ge]0,g [, so that R(g) = , R(g) = . 
We insert into (1681) and obtain 



Since R(g) ^ 0, w(g) must be equal to zero. 
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q.e.d. 

To provide an overall view, we summarize our results: 

• All regular solutions of w lie between +1 and — 1. 

• Nontrivial, regular solutions exist only for A / 0. 

• If we exclude negative values of R, it has only one maximum. 

• R can have saddle points; at these points, w has an extremum. 

All these results are of particular interest when we know whether the solutions 
have any symmetry. For example, a reflection at the equator would be nice. 
We will now consider the reflection at an axis and the point reflection (G = 
1). g then runs from to g . Because of the symmetry the interval is 
divided into two new ones. We don't want to assess where this is, but chose 
arbitrarily: 

g E [0, g ) =: I 

£':=f, Q"-=Q -^e, eeR + \{0} • (87) 

This leads us to 

: h 



g'E [0,fj =:h, q" £ 
h U h = I ■ 







Our symmetries shall be: 

1) >,C ) i/\ ( (reflection) (88) 

w{g ) = w[g ) I v ' y ' 

R(g') = R(g") (reflection) 1 
w(g') = —w(g") (point reflection) J 

First, we chose e — 2. Thus, symmetry 2) describes the properties of the 
solution of Ding & Hosoya f2~8f , equations (|7^) , (|75|) . 
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From (P 


7|) and (| 


|) we have 










Rio') 


= Rip") 


(90) 






R(e') 




(91) 






R(g') 


= W') 


(92) 






lu(g') 


= ~w{q") 


(93) 








= wig") 


(94) 








= -w(g") 


(95) 






vV) 


= V(g") 


(96) 








= -\v iw {q") 


(97) 






70') 


= -lid') 


(98) 








= i{e")- 


(99) 



On the other hand, the conditions ©, (0)' (0)' ©> (EZD > ©> 

upon insertion into ( |67|) - flBTf ), must lead to the equations (|9~2|) , ([95]) , (|99f) . 
We will execute: 



w(^) 



^K^(^) + \w{Q')l{Q.') = 

~v iw {g") - ^w(q")i(q") = -w{e") 

--J2( e ) 7 ( e )-G_ = 
-I%")7(^)-G^ = 



2 1{g) R 2 (g') 



2GV(q') + R 2 {q') 



R 2 (g"^ 



2GV{g") + R 2 {g"" 



7(0") ■ 



Symmetry 1) is also easily explained. It corresponds to the solution I have 
found numerically, as we will see in chapter [|. 

Now, we take an arbitrary e. As before, we obtain: 



R(g') = R{q") 



(100) 
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R(g') = 


-(s - l)R(n") 


R(n') - 


( F _ l) 2 fi(n") 


w(g') = 


-w(g") 


w(g') = 


(e - l)w{g") 


w(g') = 


-(£ - l) 2 w(g") 


V(g') = 


V(g") 




-^V lW (g") 


i{d) = 


1 


(e - l)R(g")R(g") 



(101) 
(102) 
(103) 
(104) 
(105) 
(106) 

(107) 



2G(V(g")-(e-l) 2 w 2 (g")) + 



+ ({e-l) 2 R 2 {g")-l)+AR 2 {g" 



(108) 



If we insert (|TD§, ([107D , (|T08D into (|7]), we see that an identity with QTUg ) 
can be achieved only if e = 2 applies. We don't need to calculate any further. 
It won't help if flB8| ) and ( j69| ) apply for an arbitrary e, since (|6"7|) doesn't. 
And we have already observed the case e = 2. 
Summary: 



R(g - 



w 



(|) = - W (g Q - |) «,(*) 



i?(^ - 
U7(^ - 



G [0, . 



Together with Prop. § and Prop. ||, we are now able to sketch all possible 
solutions. 
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3 Numerics 



As we tried to solve the system (|67|) - (|70f) by means of a standard routine for 
systems of ordinary differential equations, we concluded that most routines 
were unfit to do so. Only two of them were satisfactory: D02BAF from the 
NAG library and a routine taken from Ref. ||33|| . 

I tried to solve the equations (|6T| ) - ( |73|) with a standard shooting pro- 
cedure for solving two-point boundary value problems J33J, using both [33 



and D02BAF. In the program taken from [Q, I just replaced the routine 
for solving differential equations, leaving the rest of the program unchanged. 
The idea of the method mentioned in [R3J] is to find such initial values ( |7Tf ) - 
fl82"|) that the solutions meet at any point between the two boundary values. 
Unfortunately, it turned out that the procedure was only satisfactory enough 
when the initial values (fT9|) - fl32|) were close enough to the solution. At a 
distance of 10 -2 from the actual values, the computer was no longer able to 
find the solution. I tried this with the already known solution fl74|), ([751). 

So I returned to the standard program. I constructed the initial values 
with the expansions ( ]75| ) - (Pop, at g = 1CT 5 . By varying the parameters b 
and A, I shot for global solutions. The tolerance was of 10~ 8 . The program 
was supposed to run until R = or \w\ > 1. Now, I proceded as follows: I 
first shot from the North Pole (i.e. at g = 0) and directed b and A until the 
solution became regular. The initial values came from the conditions (|79|), 
fl80|), (|83|), (jS5|). The program told me now at which g the South Pole was, 
i.e. how big g was. Then I reversed the routine and shot from the South 
Pole. The necessary initial values now came from (PT|), (]S3), ( [SID and (PS)), b 
and A were virtually identical. It is not until we obtain the same solutions by 
shooting from both sides that we can be sure to have found a real solution. 
At this point I also employed the routine taken from |34| , as the values were 
close enough to the actual solution. For this program, it was necessary to 
introduce a new independent variable t, so that g became a function. I 
defined: 



g( t ) := gt, tE [0,1], 

e(o) = o, 

q(1) = g . 

Thus, now there are six parameters; namely b, A, g both for the North and 
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the South Pole. Now I approximated very accurately the parameters to the 
actual values. 
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4 Results 



arc 



As we have seen, the parameters of the solution of Ding & Hosoya |28 

b = 0.25 
A = 0.75 

g Q = y/2 7i , 

for G = 2. The functions R, w, 7, L are shown in Fig. [I]. 

If we assume A = 0.75 and increase b, w does not reach the value —1; it 
rather moves towards zero again. An extremum is formed. If we decrease b, 
w becomes smaller than —1 and disappears, according to Prop. [I], to infinity. 

If we assume a A other than 0.75 and try to arrange b so as to allow a 
solution with a zero in w, we can observe the following two cases: 

A > 0.75: In the best-case scenario, we can arrange b such as to make R and 
w look rather good. On the South Pole, however, there is a peak that 
grows with an increasing A. In Fig. pj, you can recognize the peak. In 
the graph, R has been drawn against R. 

A < 0.75: In this case, there is a indentation at the South Pole. Again, I 
have plotted R against R (see Fig. ||). 

In both cases the following applies: the farther away A is from 0.75 the 
farther the abscissa of the point with the ordinate —1 moves off from zero. 
Thus, I assert that for a zero of w, there is only one A. In other words: for 
one zero, there is only one regular solution. Now the question arises whether 
there are also regular solutions with more than one zero in w, and whether 
there is also just one regular solution for any number of zeros. From the 
analogy with the asymptotically flat case [|J, we can expect that. And in 
fact, the same behavior was shown for two zeros, as described above. I found 
a solution with the parameters 

b = 0.42976 
A = 0.364 
g = 8.64. 

The solution is shown in Fig. f|. We see that the symmetry 1) is fulfilled. 
There seem to be no regular solutions with two zeros for other values of A. 
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0.3 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

r/Vc r/Vg 



Figure 2: R versus R with b = 0.23655, A = 0.8 
Figure 3: R versus R with b = 0.2640625, A = 0.7 



So far, I could not find any solutions with more than 2 zeros, either. We 
presume that either the numerics fail or a function in the system becomes 
singular. As an example, we observe that the functions 



w 
R 

7 
L 



1 




cos 



solve our system (|67| ) - ([701). 7 becomes, however, singular at the equator. 
Thus, in our coordinates we don't find the solution on the computer, although 
it solves our system. That could mean that we cannot find all solutions in 
these coordinates. We must try to introduce more appropriate coordinates. 
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Figure 4: The solution with b = 0.42976, A = 0.364, g = 8.64 
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5 Time dependent cosmological solutions to 
the EYM system 

Our metric shall be 



9 = a\t) 



-dt 2 +h 



(109) 



with h being the standard metric on S 3 . The radius will be taken equal to 
1, which is always possible. We choose the gauge potential 



A = f9 



(110) 



where 9 is the Maurer-Cartan form on S 3 , understood as SU(2). / shall 
depend exclusively upon time: / = f(t). 

We set up the coupled field equations: 
For the YM equations we can use the metric g = (—dt 2 + h), since they are 
conformally invariant. It is helpful to express the metric h in terms of 9. Let 



9 
9° 

then h 



dt ; 



(r = a/i) 



With the help of the Maurer-Cartan equations 

d6 + 9 A 9 = d9 
we find the field strenght tensor 



1, 



or 



F = f9°A9 + (f 2 -f)-[9,9] 

f9°A9 1 r l + f9°A9 2 r 2 + f9°A9 3 r 3 
+2{f 2 - f) [0 2 A0 3 n + 9 3 A9 1 t 2 + 9 1 A9 2 t 3 



(111) 

(112) 
(113) 

(114) 
(115) 

(116) 



The upper line of (|1 16|) conteins the electrical components, and the lower one 
the magnetic components. Now, ( 116 ) gives us very quickly 



(F\F) 



-/ 2 + 4(/ 2 -/ 



•3 



(117) 
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Figure 5: The potential V(f) 
Figure 6: The phase portrait with E = 0.1, 0.2, 0.25, 0.3 



Here we have used the following normalization of the scalar product for the 
Lie algebra 

(X,Y) = - l -tr{X-Y) (i.e. <t;,7}> = <%) . 

Thus, the YM equations are reduced to a 1-dimensional "mechanical pro- 
blem" with the Lagrange function 

L YM = f 2 -4(f 2 -ff = T-V. (118) 
The corresponding energy is T + V = E: 

f 2 +A(f-f) 2 = E. (119) 

In Fig. [|, we plot the potential V and also the phase portrait (Fig. ||). There 
are two areas of stability at / = and / = 1, while exactly in the middle, at 
/ = 1/2, we have an unstable point, which turn out to be gauge equivalent to 
the DH solution as already mentioned in the introduction. Now, we explain 
this in short. We define the following map 

g:S 3 ^SU(2), x^x 4 -l + ixa, 

iGR 4 , \x\ = 1 , 
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and the one-form 

= g~ l dg ■ 

We choose the adequate coordinates 



r = \x\ , x 4 = y/l — r 2 . 
Then the following applies 

9 = (x A • 1 — % a x^j (dx A + ia dx 
After a short calculation we obtain 



— 9 — & j Vl — r 2 dx + 



X 



x A dx 



and finally (r = a ji) 

—9 = — ^=^= r r dr + r y/l — r 2 (t$ dd + t v sin-# dip) 
V 1 — r 2 v ' 

—r 2 (r^ dd — r$ sin-# dipj . 
If we define r = sinx, it follows that 

—9 = r r dx + sinx cosx (tq dfi + r v sin# dip^j 
— sin 2 x (t<p dd — t# sin-# dip^j . 

The BK ansatz is 



A = 



w — 1 



eh? — t,? sin-# d<£> 



We transform this with G = cosx + r r si n X an d obtain (A = 2%): 

G~ l AG + CT 1 = - r r rfA + ^ C ° SA - 1 (r y rffl - r tf sim? rfy>) 

+ — sin A (r.# di? + t v sin-# d</?) . 

This is equal to \ 9 for w = cosA and x = ^ ■ Thus, we find indeed 
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We set up the Einstein equations with A. Since tr(T) = 0, we consider 
first the trace of the Einstein equations: R = 4A. Since 



m = 4 



a -+i 

a 



we obtain 



2A n 

d + a = — a 
3 



(120) 



With respect to the orthonormalized tetrade {a0} of g, G 00 = ■% 
The energy density is given by (use ( |116| )): 



87rT f 



00 



f 2 



4(/ 2 -/) 2 



• 3 



-AT + V). 



(121) 



Thus, the Friedman equation with the A-term is: 

a* + a 2 = T + V + -a 4 . 

3 

For every T + V = E, this is again a "mechanical problem" for a(t): 

A 



(122) 



(123) 



a 2 +U(a) = E, U(a)=a 2 a 4 . 

3 

Let us first consider the case A = 0. Then 

a + a = , d 2 + a 2 = E , 

hence 

a(t) = a smt. (124) 
When we take / as being static, the derivative of ( |119[ ) will lead us to: 

/(/-1)(2/-1) = 0. 

We obtain three solutions: 

1 



/x = 0, f 2 = l, h 



(125) 
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Figure 7: The potential U(a) with A = 3/4 



For /i, A equals zero, and /2 corresponds to a pure gauge A = 9. For the 
interesting case / = 1/2, _E equals 1/4. If we shift to the physical time t, 



There is a static solution for A ^ 0, namely / = 1/2, a = (2 A/3) -1 / 2 (this 
follows from ( |120[ )). In Fig. |7|, we see the shape of U(a). The a above is 
exactly the critical point of 17(a). The DH solution corresponds to the local 
maximum of the potential V(f) and is therefore unstable. 

A stability analysis of the new solution, as well as other material will be 
published elsewhere f35fl . 
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